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ABSTRACT 

We present a critical analysis of the observational constraints on, and of the theoretical modeling 
of, aspects of cosmic ray (CR) generation and propagation in the Galaxy, which are relevant for the 
interpretation of recent positron and anti-proton measurements. We give simple, analytic, model 
independent expressions for the secondary p flux, and an upper limit for the secondary e + flux, 
obtained by neglecting e + radiative losses, e + / (e + +e _ ) < 0.2±0.1 up to ~ 300 GeV. These expressions 
are completely determined by the rigidity dependent grammage, which is measured from stable CR 
secondaries up to ~ 150 GeV/nuc, and by nuclear cross sections measured in the laboratory, p and e + 
measurements, available up to ~ 100 GeV, are consistent with these estimates, implying that there is 
no need for new, non-secondary, p or e + sources. The radiative loss suppression factor f Sif ,+ of the e + 
flux depends on the e + propagation in the Galaxy, which is not understood theoretically. A rough, 
model independent estimate of f se + ~ 1/3 can be obtained at a single energy, e ~ 20 GeV, from 
unstable secondary decay and is found to be consistent with e + measurements, including the positron 
fraction measured by PAMELA. We show that specific detailed models, that agree with compositional 
CR data, agree with our simple expressions for the e + and p flux, and that the claims that the positron 
fraction measured by PAMELA requires new primary e + sources are based on assumptions, that are 
not supported by observations. If PAMELA results are correct, they suggest that f s e +(e) is slightly 
increasing with energy, which provides an interesting constraint on CR propagation models. We argue 
that measurements of the e + to p ratio are more useful for challenging secondary production models 
than the e + /(e + + e~) fraction. 

Subject headings: Astrophysics - High Energy Astrophysical Phenomena, High Energy Physics - Phe- 
nomenology 
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1. INTRODUCTION 



Recen tly, the positron dAdriani et~al . 2009a) and anti- 
proton (Adri ani et al.l l2009bD Cosmic Rav (CR) fluxes 
were measured to high energies, e ~ 100 GeV, with un- 
precedented accuracy by the PAMELA satellite-borne 
experiment. An anomalous positron overabundance, 
compared to the expected abundance of secondary 
positrons produced by interactions of CRs and ambi- 
ent (inter-stellar) nuclei, was reported to exist in the 
energy range 1.5 GeV < e < 100 GeV based on the re- 
ported ratio of the e + flux t o the sum of the e~ and e + 
fluxes (j Adriani et al.l [200 9a). This reported overabun- 
dance was widely claimed to neces sarily imply the exis- 
tence of primary e + sources (e.g. lAdriani et al.l r2009a; 
iMorselli fe Moskalenkcj I2008D . Both astrophysical ob- 
iects (see e.g. IProfumoll2008l and references within) and 
ann ihilation of dark m atter particles in our Galaxy (see 
e.g. I Meade et~a l. 2009, and references within) were sug- 
gested as possible sources. 

Assuming that the PAMELA m easurements of the 
positron fraction are correct (see e.g. iSchubnell 2 0091 for 
cautionary notes) , the robustness of the claim of the ne- 
cessity of a primary e + source depends on the robustness 
of the theoretical predictions for the secondary e + flux 
and on the reliability of the estimate of the e + -I- e~ flux 
theoretical or observational. As for the e + + e~ flux 
recent precise measurem ents are available fflAbdo et al 
120091; IChang et al.ll2008D and also ([Aharonian et alJ2008t 
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IH. E. S. S. Collaboration: F. AharonTanl [2009]) ], espe- 
cially in the high energy range e > 30 GeV. As for the 
theoretical predictions of the e + flux, these were obtained 
in the majority of recent publications by using detailed, 
often numerical calculations of specific models, that are 
based on differing assumpti ons regarding CR generation 
and propagation (se e e.g. IMorselli fc Moskal enkol 120081 : 
IDelahave et al.l [20081 and references within). The most 
commonly used, unestablished assumption is that the CR 
propagation is described on large scales by an isotropic 
diffusion equation. Moreover, in most cases the diffusion 
coefficient is assumed to be homogenous or to exhibit a 
simple distribution with energy independent boundary 
conditions. In these models, the known relevant propa- 
gation properties are mixed with unestablished assump- 
tions and it is difficult to separate between robust and 
model dependent results. 

In fact, a rather accurate estimate of the flux of any 
secondary CR resulting from CR-Inter Stellar Medium 
(ISM) interactions, can be made wi thout a detailed un- 
derst anding of the CR propagation (Gai sser fc Schaeferl 
1992). The spatial and temporal dependence of the 
source functions for all such secondaries are the same, 
as they are proportional to the product of the CR in- 
tensity and the ISM density (assuming that the primary 
CRs have a uniform composition) . The ratio of the local 
fluxes of two secondary particles should be equal to the 
measured ratio of their local production rates. Thus, a 
measurement of one secondary specie (e.g. Boron), al- 
lows a prediction of all other secondaries, given that the 
interaction cross sections are known. The main challenge 
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in the analysis comes from the fact that CRs suffer sig- 
nificant losses, which are different for different species, 
during their propagation. Nuclei suffer spallation and 
decay, while positrons suffer energy losses due to Inverse 
Compton and Synchrotron emission. Anti-protons suffer 
some losses due to annihilation and scattering. 

Guided by these considerations, we present in this pa- 
per a critical analysis of the observational constraints on, 
and the theoretical modeling of, aspects of CR genera- 
tion and propagation in the Galaxy, which are relevant 
for the interpretation of recent e + and p measurements. 
The plan of the paper, along with some of the principle 
points, is as follows. 

In § [2] we review the main observationally established 
properties of CR propagation. We focus on the analysis 
of stable CR nuclei measurements and show that they 
can be accurately modeled using transparent physical ar- 
guments. In particular, the losses due to spallation are 
addressed. The CR grammage, extracted from the stable 
CR nuclei measurements, allows us to write down in § [3] 
analytic, model independent constraints for the expected 
secondary p and e + fluxes, including an accurate estimate 
of the p flux and an upper limit for the e + flux, obtained 
by neglecting the e + radiative losses. Anti-proton and e + 
measurements available up to e ~ 100 GeV, including the 
PAMELA positron fraction measurement, are shown to 
be consistent with these estimates. 

The e + flux measurement can be used to infer the sup- 
pression factor / Sj e+(e) due to radiative energy losses, 
defined by the relation 



(^) /s,e+ (^)^e+ , no losses (^) i 



(1) 



where n e + no i oascs is the denisty of positrons that would 
have been observed if positrons did not suffer energy 
losses. Similarly to secondary nuclei and anti-protons, 
n e+,no losses can be calculated directly using the measured 
CR grammage. We show that the existing e~ and e + 
data imply that f s , e +(s) ~ 1/few in the PAMELA en- 
ergy range, and is indicated to be slightly growing with 
energy. 

In §|4]we discuss the measurements that allow an indi- 
rect estimate of f s e +(e), namely via comparison to the 
suppression of the abundance of radioactive elements. 
We show that observations allow a rough, model in- 
dependent, estimate of f Sie + ~ 1/3 at a single energy 
e ~ 20 GeV. This estimate is shown to be consistent 
with e + measurements, implying that (i) positrons are 
almost certainly entirely secondary, and (ii) e + measure- 
ments at high energies e > 20 GeV provide novel infor- 
mation on the temporal behavior of CR propagation that 
is not currently accessible by other types of experiments. 

Any propagation model that reproduces the CR com- 
positional measurements discussed above will be sub- 
ject to the constraints we provided for the positron and 
anti-proton fluxes. In particular it should agree with 
our estimates for the secondary anti-proton flux, should 
agree with the estimate for the secondary positron flux 
at e ~ 20 GeV and should give positron fluxes lower 
than the upper limit we derived, at positron energies 
e > 20 GeV. In section § [5] we illustrate this point by 
briefly discussing specific, commonly used models for CR 
propagation, focusing on the Leaky Box and disc+halo 
diffusion models. Details are given in appendix § [AJ 



The equations for stable nuclei and anti-protons in the 
Leaky Box model and thin-disc+halo diffusion models, 
are shown to be equivalent to our equations, as long as 
the grammage is set to be equal to the measured gram- 
mage. We reproduce the results of detailed calculations 
of anti-protons and positrons using simple analytical ap- 
proximations for the grammage and energy loss suppres- 
sion factor. We note that similar simple analytic ex- 
pressions for the expected fluxes of secondary particles 
in th ese specific models can be found in the literature 
(e.g. iGinzburg fe Ptuskinl 119761: IBerezinskii et al.l 119901 : 
ILongair. M. SJI1992t IPtuskin et all 1200911 . We highlight 
the main, frequently unestablished, assumptions leading 
to the claimed conclusion that secondary models fail to 
reproduce the e + measurements. 

Previous claims, that a primary e + source is necessary 
to explain the positron fraction observations, are dis- 
cussed in § [6] We show that these claims are based on as- 
sumptions regarding the generation spectrum of primary 
electrons and regarding the propagation of electrons and 
positrons, which are not supported by observations. 

In section §[7] we discuss specific models, in which the 
positrons are of secondary origin, that were s uggested to 
explain recent e + abundance measurements (|Blasill2009l : 
IShaviv et al.l [20091: ICowsik fc Burehl [2009ah . Although 
these models do not adopt the assumptions discussed in 
§ HI which lead to a conflict with the PAMELA results, 
they do adopt other assumptions, which lead to conflicts 
with constraints derived from other CR measurements. 
We emphasize, though, that there is a variety of mod- 
els that can be constructed, in which both e + and anti- 
protons are of secondary origin, and which are consistent 
with the constraints arising from CR observations. The 
main point that should be addressed in such models is 
the energy dependence of the e + suppression factor f SyC + 
at energies exceeding e > 20 GeV, which is indicated by 
the PAMELA measurements to be slightly rising with 
energy. 

Our results are discussed and summarized in § [8] 
Throughout the paper, we limit the discussion to rel- 
ativistic e > 10 GeV/nuc energies, since at sub GeV 
energies the theory of CRs becomes more complicated 
and uncertain. This is due to the presence of various ef- 
fects, including ionization losses, energy dependent spal- 
lation cross sections and charge dependent solar modu- 
lation etc. Some of the observed effects are not well un- 
derstood (e.g. solar modulation) and possible additional 
effects may exist (e.g. reacceleration). 

2. CR GRAMMAGE 

A complete coverage of the theoretical and observa- 
tional efforts to understand CR origin and propagation 
is beyond the scope of this paper. In this section and in 
section § |4] we emphasize what we believe are the main, 
observationally confirmed properties of CRs that allow 
for predictions of secondary e + and p fluxes. 

Although the acceleration and propagation of CRs 
in our Galaxy has been studied for many decades 
Ginzburg fc Svrovatskiil 
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eil | (|1987t ): IBerezinskii et al.l ... . 
)92|) and more recently iStrong et all 

to 

date still only little is known. Current knowledge is 
essentially restricted to the following features: 



Blandford & Eichle: 
Lon gair. M. S.I (|199 
420071) ]. Dased on results of numerous experiments 
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1. A qualitative global picture. CRs are generated 
in the Galaxy and are confined by the non trivial 
magnetic field of the Galaxy to long, rigidity depen- 
dent time scales (~ 10 7 yr at e ~ 1 GeV). The CR 
halo is probably considerably thicker than the few- 
hundred-pc-thick gaseous disc. Indirect evidence 
(notably synchrotron emission from nearby galax- 
ies and gamma-rays from our Galaxy) suggest that 
the CR distribution is, to an order of magnitude, 
homogenous throughout the halo. 

2. Some quantitative features of CR propagation, es- 
tablished by CR compositional measurements. The 
most important quantitative feature relevant to 
the calculation of secondary particles is the rigid- 
ity dependent average amount of column density 
(grammage) traversed by CRs. Information on CR 
grammage is based on measurements of the rela- 
tive abundance of spallation generated secondary 
particles such as Li, Be and B or Sc,Ti and V up 
to energies e < 150 GeV/nuc. Limited constraints 
on the propagation time scales exist due to the ob- 
served suppression in the abundance of unstable 
radioactive nuclei, especially 10 Be at low energies, 
e ~ 100 MeV/nuc. 

As CRs traverse the interstellar medium they suffer 
spallation. This leads to a few features that are observed 
in the local CR nuclei spectrum and composition: 

• An over abundance with respect to the solar com- 
position of chemical elements that are produced in 
the spallation interactions. For some elements like 
Li,Be,B and Sc,V,Ti, this over abundance is by a 
few orders of magnitude, and thus these are most 
probably purely secondary particles. 

• A suppression of the flux of particles that were lost 
due to spallation. This effect grows with the total 
cross section for spallation, which in turn grows 
with the particles mass and reaches a factor of a 
few for Iron. 

The "optical depth" for spallation of a particle is pro- 
portional to the amount of grammage surpassed by the 
particle, irrespective of the time it took the particle to 
accumulate this grammage or to the densities through 
which it traversed. 

In this section, we discuss the measured amount of 
grammage deduced from the compositional analysis of 
stable nuclei. This allows us to derive in § [3] model in- 
dependent expressions for the secondary p flux, and an 
upper limit for the e + flux, obtained by neglecting e + 
radiative losses. For clarity we begin in § 12.11 bv consid- 
ering secondary particles that do not suffer losses. The 
definition of the CR grammage and the use of secondary 
particles to deduce it's value are discussed. Spallation 
losses are addressed in § 12.21 where we discuss the obser- 
vations and the measured grammage values. The analysis 
of measurements of radioactive unstable nuclei are dis- 
cussed in § HI where they are used to derive constraints 
on the e + radiative losses suppression factor. 

2.1. Stable secondaries without spallation losses 



Consider hypothetical CR secondaries that suffer no 
losses at all. The local density of such CRs at a given 
rigidity will be proportional to their local generation rate 
under the following assumptions: 

(2) 

1. secondary particles with the same rigidity propa- 
gate through the ISM in the same way (diffusively 
or otherwise), 

2. the rigidity of the products equals the rigidity of 
the primary, 

3. the energy of secondary particles does not change 
during propagation, 

4. the composition (but not necessarily flux, spectrum 
or target density) of CRs is uniform throughout the 
region in which most of the secondaries observed 
here are produced and during the time they were 
produced. 

Under the above conditions, for any two stable CRs A 
and B at a given rigidity, the following equation holds: 

n A _ Q A ,„.. 
riB Qb 

where Qi is the local production rate of the secondary i 
and is given by 

Qi = nj^^cpisM, (4) 

where Cj-m is the decayed spallation cross section of the 
parent nuclei j into the secondary i per ISM nucleon. To 
see this, note that since the composition is the same at 
any given time or place, the ratio of the production rates 
will be the same everywhere. For every particle B that is 
generated, Qa/Qb particles of type A are generated. As 
the particles propagate in the same manner, the density 
of A will be Qa/Qb times that of B. 
It is useful to write equation (O as 

ni {e) = ^Mflx osc (e/Z), (5) 

PISMC 

where the grammage A osc , defined by this equation, pa- 
rameterizes the column density of target material tra- 
versed by the CRs and is the same for all species. 
Qi(e) / (pismc) is the local net generation of the CRs per 
unit traversed mass, and is independent of the local ISM 
density. Note that rij are directly measured by CR ex- 
periments and o~j , o-j^i are measured in the laboratory. 
Thus, A eS c can be directly extracted from compositional 
measurements of CRs using Eqs. (HI and Once X esc 
is determined, the density of any secondary can be com- 
puted using these equations along with the primary CR 
and the cross section measurements. 

2.2. Stable secondaries that suffer spallation losses 

In reality, all measured spallation secondary nuclei suf- 
fer significant losses due to spallation. Equations (JS])-© 
can be generalized to include spallation by replacing the 
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production rates Qi (e)/ (pismc) with net production rates 
Qi given by 

- Q t n l a l 

Qi = , (6) 

Pismc m p 

where Ui is the cross section for destruction of the CR 
per ISM nucleon and is approximately given by rj 
40A^ 7 nib, with Ai the atomic number of the particle i. 
This prescription results with 



n B 



(7) 



and 



«i(e) = Q z (e)X esc (e/Z). 



(8) 

To see this, note that the distribution of particles A that 
suffer spallation losses is equal to the distribution of par- 
ticles A' that do not suffer spallation if the production 
of A' is equal at any given time and place to the net 
production of A. 

Equations © and ([8]) can be written in a directly ap- 
plicable form as: 



X csc Qi/(pisMc) 



1 



(9) 



As far as we know, Eq. © is consistent with all data 
of CR composition. Above a few GeV/nuc, the value of 
A csc is measured to be 



8.7 



(lOZ Gev) 



-0.5 



g cm 



(10) 



with different fits varying by ~ 30 % in the range 
10 GeV < e/Z < 100 GeV (e g lEngelman et~al~l fisM 
Uones et al.ll200lUWebb"er~ll2003l ). There are indications 
that the power law behavior of X Rsr continues to hun- 
dreds of GeV (IBinns et al.|[T98l lAlm et al.ll200l . see 
however ([Zatsepin et al.ll2009f ). Henceforth, we assume 
that the grammage parametrization given in Eq. (|10|) 
holds up to e/Z ~ 300 GeV. 

For strong spallation losses, a 3> m p /X esc , the density 
of secondaries given by Eq. © approaches a value that 
is independent of A osc , 



p 



(11) 



The strong suppression due to spallation of the heavy 
secondaries (Sc, V, Ti), results in a small deviation of 
their densities from the limit of infinite grammage 



X e 



0.3 



-0.7 



10 GcVZ 



0.5 



Thus, the measurements of these sub-iron elements are 
useful for determining A csc only at high energies e/Z > 
100 GeV. 

What makes equation 0-© non trivial, is the fact 
that the loss term has to be included in the expression 
for the net generation rate, Eq. © . As the net generation 
rate of two particle species A and B is affected by their 
own density, their relative abundance is required to be 
uniform in order that Qa/Qb be uniform. Thus the 



validity of equations ©-© suggests that the relative 
abundance of the secondaries themselves is uniform. 

Perhaps, the simplest propagation model in which the 
above conditions are realized is the homogenous Leaky 
Box Model (LBM, see §E]). Obviously this model sat- 
isfies the conditions © and thus equations ©-© are 
guaranteed to hold. It is also known that these equa- 
tions are satisfied for disc-halo diffusion models in which 
the radial extent of the cosmic ray halo is much larger 
than the scale height, which i s in turn much larger than 
the width of the gas disc (e .g. iGinzburg fc Ptuskinlll976t 
Schlickcis er fc Lerchdll985fl . In this case again, condi- 
tions © are trivially satisfied. Expressions for A csc in 
these models are given in § [5] and § O 

In our view, Equations © - © are natural relations, 
that are expected for a wide range of models, that sat- 
isfy conditions © and that were empirically validated. 
Diffusion models, that have a large halo and thin disc 
(virtually all currently used models) or leaky box models 
are particular models that satisfy these conditions. In 
fact, these equations will hold in any 1-D model in which 
the gas is concentrated in a thin disc, assuming that the 
transport of particles (diffusive or otherwise) depends on 
their rigidity only and that their energies do not change. 

The good news are that these equations allow us to ob- 
tain robust predictions for secondary particles like anti- 
protons and positrons (with the latter requiring more 
care due to energy losses, see § [3] and § 0J. The bad 
news are that the measurements of stable secondaries, 
which are probably the single most important type of 
measurements for quantitative research of CR propaga- 
tion, carry little information regarding the precise form 
in which CRs propagate. 

We conclude this section with a comment about the ap- 
plication of these equations to primary CRs. Unlike sec- 
ondaries, the source function of primary particles is not 
known. In fact, equation © is used to deduce the aver- 
aged source spectrum and in particular the total required 
energy output of CRs. The fact that the resulting source 
spectrum, when using the grammage deduced from the 
secondary measurements, i s approximately the sam e for 
the different elements (e.g. Engelma n et al. Ill990h sug- 
gests that the same equations are applicable to the pri- 
mary CRs as well. This in turn suggests that the prop- 
agation of CRs averages the generation spectrum over 
distances larger than the inhomogencitics of the primary 
sources. We note that the application of these equations 
to primaries is somewhat less substantiated theoretically 
and observationally than for the secondaries. 

3. APPLICATION TO SECONDARY 
ANTI-PROTONS AND POSITRONS 

In this section, we estimate the expected flux of 
positrons and anti-protons using the measured CR tra- 
versed grammage discussed in § O We first discuss the 
local production rates of positrons and anti-protons in 
§ 13.11 We then write down a model independent expected 
p flux in § 13.21 and show that it agrees with observations. 
An upper limit for the secondary e + flux, obtained by 
neglecting energy losses, is given in § 13.31 and compared 
to observations. The e + energy losses are addressed in 

§m 

We note that for positrons and anti-protons the sec- 
ond of conditions © is not satisfied, as these particles 
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are generated at rigidities lower by a factor of ~ 10 com- 
pared to their progenitors. The fluxes of positrons and 10° 
anti-protons are therefore more sensitive to spectral vari- 
ations of the CRs in the Galaxy compared to products 
of spallation of nuclei. 



3.1. Production of anti-protons and positrons 

The rate of production of positrons and anti-protons 
depends on the flux of primary CRs, the ISM composi- 
tion and nuclear cross section data. For e > few GeV, the 
cross section dependence on energy is essentially dictated 
by the cross section for pp collisions, while the presence of 
heavy target and projectile nuclei can be approximated 
via an energy independent scaling factor, denoted here 
by £s,A>i, of order unity. 

Concerning the production of ant i-protons, we adop t 
the cross section parametrization of iTan fe Ng (1983b). 
Given a measurement of the primary proton flux J p , the 
p production rate per unit ISM particle mass is given by 



poo 

Qp(e) = 2£p ij4>:L 47r / de p J p (e p ) 

J e-n 



da p (e p ,e) 
de v 



(13) 



where the factor of 2 accounts for the decay of antineu- 
trons produced in the same interactions. 

For positrons, the production rate is given by a for- 
mula similar to (|13p. using the cross section for fi- 
nal state positrons resulting from the decay of charged 
mesons. For the charged mes on cross section, w e again 
adopt the parametrization of ITan fc Ngl (1983b). The 
subsequent e + yield is calculated using standard elec- 
troweak t heory. The e + yield we find agrees with the 
results of iDelahave et al.l $2008) to - 10% for the same 
parametrization. 

For a steeply declining, smooth primary proton spec- 
trum, it is useful to parameterize the resulting local gen- 
eration rate of anti-protons and of positrons per ISM 
mass by the following equation: 



sQs(s) = ts,A>i(e)C s , pp (e)4Tr(10£)J p (l0e) 



0~ppfl 



(14) 



where S = p,e + stands for anti-protons and positrons 
respectively, <7 PPj o = 30 mb is a cross section normaliza- 
tion chosen to be approximately the inelastic cross sec- 
tion for pp in teractions at the energy range 10 GeV < 
e < 300 GeV (|Tan fe Ndll983al) and C s , PP {e) is a dimen- 
sionless coefficient that weakly depends on the primary 
spectrum. 

The values of C PtPP and C e + pp are shown in figure [T] 
for a power law proton flux J p oc e -7 with 2.6 < 7 < 2.9. 
For this range of values for 7, the approximation 



C PtPP = 0.18 - 0.041ogi (e/500 GeV), (15) 



and 



C e +, PP = 0.58 



(16) 



at the p (e + ) kinetic energies 10 GeV < e < TeV, is 
accurate to better than 10% (approximation shown in 
the figure). 

For completeness, we also give the approximate 
value of C e + i)9 p corresponding to the parametrization of 




Kinetic Energy [GeV] 

Fig. 1. — Production coefficients Cs, P p(s) [Eq. JH}] of anti- 
protons and positrons. The lower blue curves are the production 
coefficients of anti-protons for proton spectra 7 = 2.6 — 2.9, and 
the upper red curves are the production coefficients of positrons for 
the same valu es o f 7. The black dotted line is the approximation 
given in Eqs. II15I I and dlTI ■ 



IKamae et all (|2006[ ) 

C e + pp = 0.31 + 0.15 log 10 (e/100 GeV) 9 (e - 100 GeV), 

(17) 

valid in the range 2.6 < 7 < 2.9 and 10 GeV < e < 
1 TeV to an accuracy of 10%. 

3.2. Anti-Protons 

Using equations © and (|14p . and assuming a power- 
law proton spectrum J p oc e -7 , the expected ratio of 
anti-protons to protons is given by 



■In 



-y- —10 7+1 '£p,A>iCp,p P (e 



J, 



pp,inel,0 



X, 



1 



,sc i + ^x r 



(18) 

Using the mea sured value of X P , sr , from E q . (fTOl) , s etting 
6).4> i = 1-2 (jGaisser &: Schaeferl Il992t iSimon et al.1 
fl998h . we find 

J 4- «3.6 x HT 4 ^' w(£) 



x 1 + 0.16 



0.1 j \ 

l_Op_ 

V30 m 



10 GeV 



-0.5 



mb/ VlO GeV 



-0.5 



.(19) 



In our calculation we adopt Op from lTan fc N g (1983a), 
where Op ss 30 mb holds to an accuracy of about 20% 
in the range 10 GeV < e < 100 GeV. This result is 
compared to experiments in figure [2J As can be seen, the 
results are in good ag reement with recent measurements 
(jAdriani et al.ll2009bt l lBeach et al.ll200l . 

Similar calculations, up to differences in the cross 
sect ion and grammage parameterizations, were made 
in (jGaisser fc Schaefer I 1199% ISimon et all 11998ft . In 
figure [H the results of two detailed diff usion models 
(|Donato et al.ll2009t iMoskalenko et"aL1l2002D are plotted. 
The model of ()Donato et al.l 120091 ) assumes a primary 
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secondary, ±40% band 

secondary 

" Moskalenko 2002 

• Donato 2009 



O PAMELA 
X HEAT 




1 1 1 2 

Kinnetic Energy [GeV] 



Fig. 2. — Anti-proton to proton flux ratio. The blue and dark 
brown error-bars are the Jp/J p ratios measured b y PAMELA 
(|Adriani et al.ll2009bl) and HEAT IjBeach et alj |200T). The black 
line is the expected ratio us ing equation l|18[l. with the cyan col- 
ored region denoting a 40% (Simon ct al. 1998) uncertainty band. 
The brown band d epicts the results of the diffusion model of 
IDonato et all (I2009T) . 
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Fig. 3. — I nterst ellar e+ flux measured by HEAT 
HDuVernois et al.l l200ll) c omp ared with the theoretical upper 
bound, obtained from Eq.lf23} with / + = 1 for e < 300 GeV. 
The allowed value of the e+ flux, corresponding to / e+ < 1 is 

shown as a shaded region. Also plotted i s the e+ + e~ flux as 
meas ured by HEAT IDuVernois et aL|[200ll 1. FERMI llAbdo et al.l 
[2009T) . ATIC HChang et al.H2008l) and AMSO HAlcaraz et aI.II2000l l. 
Fluxes are demodulated using the indicated values of the solar 
modulation potential <E>. 



proton sp ectral index of 7 = 2.84, cross sections ate 
taken from lBringmann fc Sal ati (2007) with the p-p cross 
sections based on lTan fc Nd (|1983aj) (which are slightly 
smaller than the parametrization by I Tan fc Ng| ["l983b, 
that we use) , and uses an energy depe ndent diffusion co- 
efficie nt D(e) oc £ 7 . The model of (jMoskalenko et all 
l2002f) assumes a primary proto n spectral i ndex o f 7 « 
2.75, pp cross sections based on lTan fc Nd ()1983b0 . and 
uses an energy dependent diffusion coefficient D(e) oc 
e - 6 . The results of the different computations agree to 
within a factor of ~ 2. 

3.3. Positrons 

Using equations CD,® and ([14| and assuming a 
power-law proton spectrum J p oc e~ J , the expected ratio 
of positrons to protons is given by 

Je+ _ f H f r .\ a PP>inel,0 v 

jp Trip 

(20) 

Using the m easured value of X esc from Eq.([T0 |) . setting 
L+,A>i ~ 1 (|Moskalenko fc Strondll998[ ). we find 

llxl0 -^(^M) (_£_)-. 

(21) 

Adopting a high energy Inter Stellar (IS) proton flux 
(jMoskalenko et al.H2002li . 

J p = 1.6 x 10 4 (e/ GeV)" 2 75 GeV" 1 m^sec" 1 sr" 1 , 

(22) 



we have 

GeV" 1 m-^ec" 1 sr" 1 . 

(23) 

An upper limit for the e + flux, obtained by setting 
f s . e + = 1 in equation (|23[) . is shown in figure [3] as a 
shaded region along with e + and e + + e~ flux measure- 
ments. As can be seen, the e + data are compatible with 
a secondary origin. The data implies a suppression fac- 
tor of order f s e + (e) ~ 1/ few (see discussion of e + energy 
losses in §[4|. 

The positron fraction was measured by several experi- 
ments including the recent PAMELA experiment. Since 
the e~ flux is poorly understood theoretically, we do not 
have a model independent prediction for the positron 
fraction. In order to compare our predictions for the 
e + flux to the positron fraction measurements we use 
e + + e~~ measurements. 

The claimed PAMELA excess extends from sub 
10 GeV energies with high statistics to 100 GeV with 
low statistics. The e + + e~ spectrum is measured by 
different experiments in different energy intervals within 
this range. Thus, it is useful to compare the PAMELA 
measurement with the various e + +e" measurements and 
the theoretical expectation ((23|) using a single plot. This 
is done in figure [H where we divide our e + flux upper 
bound derived from (|23p by the measured e + + e~ flux 
and compare the results to the PAMELA positron frac- 
tion data. 

For completeness, in figure [5] we reproduce the same 
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Fig. 4. — Interstellar positron fraction as measured by PAMELA 
(Adriani et al. 2009a) compared with an upper limit for the 
positron fraction. The upper limit is obtained by dividing the 
theoreti cal e+ flux upper boun d with e+ + e~ fluxes meas ured by 
HEAT lIDuVernois et al.ll200lft . FERMI IIAbdo et aJJl2009h . ATIC 
llChang et al.H2008l) andAMSOl (Alca raz et aTE oOO'). 
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Fig. 5. — As figur e [5] using the cross section parametrization of 
IjKamae et al.ll200g i. 



calculation leading to figure HI using th e e + produc- 
tion cross section parametrization from ()Kamae et al.l 
2006). This recent parametrization results with a e + 
yield that is lower by about a factor of two com pared 
to the parametrization from ()Tan fc Nd 1198 3b). used 
throughout this paper. 

As can be seen, the PAMELA and e + + e~ measure- 
ments are consistent with a secondary origin for the 
positrons throughout the measured energy range, with 
a reasonable suppression factor / + (e) ~ 1 — 4 (see dis- 
cussion of e + energy losses in § [4j) . 

As can be seen in figures [4] and OH the positron fraction 



cannot exceed e + /(e + + e ) sa 0.2 in the energy range 
100 GeV < e < 300 GeV. This is based on the measured 
e+ + e" flux by FERMI and ATIC, the measured CR 
grammage Eq. (fTU|) . limited by currently available data 
to hundreds of GeV, and the fact that f se + < 1. 

In the energy interval e > 30 GeV, there is an indica- 
tion from the measurements that f s e + (e) is rising with 
energy. Since f 3e + < 1, this trend cannot continue in- 
definitely. If f s . e + is rising with energy, a break in the 
rising pattern must occur at e < 100 GeV. If primary 
electrons are affected by the energy losses which affect 



the secondary positrons, a feature in the 



flux 



is expected. Note that at higher energies the statistical 
and systematic errors in the positron fraction are large, 
and thus the rising trend of f s , e + is not implied by the 
observations with high significance. 



3.4. Positrons/ Anti-Protons 

Both positrons and anti-protons at a given energy e 
are produced by the interactions of the same primary 
CRs (mainly protons) at roughly the same primary en- 
ergies (roughly lOe). Hence, the e + to p flux ratio is 
particularly insensitive to the compositional and spec- 
tral variations of the CRs. 

The ratio of positrons to anti-protons is given by: 



J r> 




c,. 



C u 



- , p,pp(s) 

(24) 

Both ^ e + j A>i an d £,p,A>i are not very different from 
1 . Measuremen ts suggest that ^ e + ! A>i < £,p,A>i 
(jAdler et al.ll2006f) . The p losses are not significant above 
10 GeV (about 15% at 10 GeV and droping with energy). 
Quite generally we thus have: 



J e + < C e + pp (e) 



■Jn 



Ip.pp 



(25) 



where Qs, PP is the generation rate of specie S in pp in- 
teractions. The e + to p branching ratio, Q e +,pp/Qp..ppi 
is calculated based on particle physics and accelerator 
measurements (see section § 13. ip and is shown in figure 
[H As seen in the figure, this ratio weakly depends on 
energy above a few tens of GeV and is given roughly 
by Q e + .pp/Qp.pp ~ 3. Furthermore, this ratio depends 
weakly on the primary proton spectrum. To illustrate 
this fact, we also plot in figure [5] the ratio resulting from 
an extreme value of the proton spectral index, 7 = 2. 
We find that Q e + , pp/Qp.pp remains within 50% from its 
value for 7 = 2.75. 

The e + to p flux ratio is not given in any experimental 
report that we are aware of. In order to compare the 
expected ratio, in the energy range e ~ 10 — 100 GeV, 
with the ratio inferred from the PAMELA measurements 
of e + /(e + + e~) and p/p combined with proton and e + + 
er measurements, it is useful to parameterize the ratio 
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Fig. 6. — e + to p production ratio due to pp int erac tions, which is 
also the flux ratio upper limit given by equation (125 1 . The different 
smooth, colored lines are for different values of the proton spectral 
index 7, in the range 2.6 < 7 < 2.9 . The cros sed black line is the 
expected ratio using equations (1 1 5 I t and H17I I. The dashed black 
line corresponds to an extreme value 7 = 2. The cyan colored 
band depicts 50% modulation over the central value for 7 = 2.75, 
to be regarded as a conservative uncertainty estimate. 



J e + _ „ ( Je+/(J e + + Je~) 
Jp \ 0.1 

Jpf^ 

1.6 x 10 4 GeV 1 ' 75 m- 2 s^ 1 sr" 1 

e x-0.25 



(J e + + J e -)e 3 



150 GeV 2 m- 2 s- 1 sr- 1 



Jp i Jp 

2 x 10-' 



30 GeV 



(26) 

All the factors in Eq. (|26p are measured to be close to 
unity in the relevant energy range. Hence the experimen- 
tal results are consistent with Eq. flMl) for / + ^2 — 3. 

Obviously, a direct measurement of the e + to p flux ratio 
would be much more useful. 

4. PROPAGATION TIME SCALES AND POSITRON 
ENERGY LOSSES 

The local density of secondaries that decay or suffer 
energy loss will be suppressed compared to the prediction 
of Eq.® if their decay/loss time is short compared to 
their typical residence time. Positrons lose energy by 
synchrotron and Inverse Compton (IC) processes with a 
cooling time 



2 4 

e mzc 



£(JtUt c 



(30 Gev) 



10 7 yr 



leV cm' 



(27) 



where Ut is the sum of the ambient photon and mag- 
netic field energy densities. This time scale is compa- 
rable to the CR propagation time scale at the energy 



range of interest, so some amount of cooling suppression 
of positrons is expected. The radiative loss suppression 
depends on details of the CR propagation in the Galaxy, 
which is not understood theoretically. 

Measurements of the abundance of radioactively un- 
stable CR nuclei with known decay times give valuable 
information on the residence time scale of the CRs and 
hence on the suppression of electrons and positrons due 
to cooling. In this section we use the results of such mea- 
surements to constrain the value of f se + [see Eq. (flj)]. 
In § 14.11 we explain how measurements of the suppres- 
sion factor of radioactive nuclei can be used to constrain 
the cooling suppression of electrons in a model insensi- 
tive way. In particular we address the complications that 
arise due to the spallation of the radioactive nuclei and 
to the difference in the suppression mechanism (decay vs. 
energy loss). In S M4.2l we discuss the results of abundance 
measurements and give an order of magnitude estimate 
of f B e + . This estimate is shown to be consistent with the 
results of § [3] A discussion of f s e + in common models is 
given in § 15.21 

4.1. Radioactive decay vs radiative energy loss 

As in the case of electrons and positrons, it is useful to 
separate the effect of decay and propagation for unstable 
nuclei. A complication that arises in the case of nuclei 
is the fact that in addition to decay, the nuclei are de- 
stroyed by spallation. It is useful to separate the effect 
of losses due to decay from that of spallation. This can 
be done in a model independent way by considering the 
net production of the particles instead of the source pro- 
duction, and defining the decay suppression factor f s j 

by 



J s.iQiX-esc fs 



i/(piSMc) - rii — 



X r , 



(28) 



where Eqs. ([6]) and ((5J were used. This suppression fac- 
tor is more useful for obtaining the suppression factor of 
positrons compared to the more commonly used surviv- 
ing fraction f s .i (e.g. iWebber fc Soutoull [T998h . defined 



as 



j s, i^-i, no decay; 



(29) 



were 7ii jno decay is the density of this nuclei that would 
result if it did not suffer decay (but did suffer spallation 
losses). Using Eq.©, we find that n,^ no decay is given by 



7 2 (Qi/PiSM)X 

C 

Js, ffii, no decay — Js,i Z 

1 - 



(30) 



Substituting n., from Eq. (|30|) in Eq. (|28|) and extracting 
f s .i, we find that the two factors are related by: 



S,l 

s.i 



1 H A csc (l — f s A 

m„ 



(31) 



Next we demonstrate that under a wide range of con- 
ditions, the suppression factor f s .i{e/Z) of a nucleous 
that suffers decay with a decay time r^, is similar to the 
suppression factor f s e + (e) that positrons suffer due to 
energy losses with a cooling time r c , if r c ~ Tg at the 
same rigidity. The loss rates due to decay and energy 
losses in a general transport equation (assuming that the 
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dtrii. 



decay 



(32) 



distribution is close to isotropic) are given by 

Td 

and 

9t «. e +, energy loss = 9 £ (sn e +) (33) 

respectively. For losses due to IC and synchrotron where 
e oc e 2 , we can rewrite Eq. ([33)l as 



d t n 



i'^e^ , energy loss 



d\og{e 2 n e + ) \ n e + 
d log(e) 



(34) 



Assuming that the e + spectrum is steeper than n oc e~ 2 , 
the term in parentheses is negative and of order unity. 
The energy loss term Eq. (f34|) is similar to the decay term 
Eq. (|32|) with Td replaced by r c , up to an order unity 
correction. In particular, for a spectrum n oc s 1 this 
term equals 7 — 2 and is unity for 7 = 3. 

The similarity between cooling and decaying can be 
understood as follows. During a time interval r c (e), a 
e + with energy e loses a significant part of its energy. 
Assuming that the e + spectrum is steeper than n oc e~ 2 , 
the contribution of such positrons to the lower energy 
flux is negligible and they are effectively lost. 

Since the loss terms of positrons and unstable nuclei 
with Td ~ t c are similar, we expect that their suppression 
factors will be similar. In l5.2l we show how this is realized 
within specific models. 

4.2. Measurements and constraints on f SyC + 

Perhaps the single most important unstable isotope 
which is used for propagation time scale estimates is 
10 Be. This is due to the fact that (a) it is a pure sec- 
ondary and (b) it has a lifetime r res t ~ 2.2 x 10 6 yr, 
which turns out to be of the same order of magnitude as 
the escape time. Other commonly used u nstable nuclei 
inclu de 14 C, 26 Al, 36 C1 and 54 Mn (see e.g. IStrong et al.l 
120071) . 

Direct measurements of the abundance of these iso- 
topes exist only at sub-GeV/nuc energies. In fact, mea- 
surements with high statistics are limited to a narrow 
energy range around ~ 100 MeV/nuc, where solar mod- 
ulation and other low energy complications are impor- 
tant. Nevertheless, these low energy data are often used 
to calibrate model parameters, which are later used to 
predict fluxes at all energies (e.g. iBerezinskii et al] | 1990 |: 
IMoskalenko fc Strond 119981 IStrong fe MoskalenkolboOli 
and references within.). 

An indirect measurement of the surviving fraction can 
be made by using the high energy, charge (as opposed 
to isotopic) composition of CRs. For example, by mea- 
suring the charge ratio Be/B, and comparing it to the 
expected ratio that would result if 10 Be decayed com- 
pletely or didn't decay at all, the suppression of 10 Be can 
be extracted. Using this method, the surviving fractions 
at e < 15 GeV/nuc (rigidity of e/e Z < 30GV), were 
meas ured for a few unstable isotopes (jWebber fc Soutoul 
1998). In particular, the surviving fractions of 10 Be 
at energies e s» 4(16) GeV/nuc [rigidities of e/eZ = 
10(40)GV] were measured to be / s ,io Be « 0.4(0.5) ± 0.1. 
Using Eqs. (fl0|) and (j3"Tj) . and adopting eno Bo w 200 mb, 
the suppression factors are found to be f s io Be w 0.3(0.4). 



The decay time of 10 Be in the observer frame at these two 
energies is Td ~ 0.88(2.2) x 10 7 yr. These decay times 
are similar to the cooling times of positrons at the same 
rigidities, t c w 3.1(0.78) x 10 7 yr(C/ T /eV) _1 , as given by 
Eq.(23. 

Similar measurements of 26 Al, 36 C1 and 54 Mn, that 
were derived based on ele mental abundances of Al CI and 
Mn, give similar results (|Webber fc Soutou]| IT998). The 
analysis using these isotopes suffers additional uncertain- 
ties since Al, CI and Mn are not purely secondary. 

Note that these measurements are particularly sensi- 
tive to uncertainties in the spallation cross sections. For 
example, the Be/B ratio changes by about 20% only be- 
tween no 10 Be decay vs complete decay. An accuracy 
of the cross sections considerably better than 10% is re- 
quired. The resulting systematic uncertainties are hard 
to estimate. The fact that the different measurements 
seem to be consistent with each other give some support 
to the derived values of f Sie +. 

Assuming that these measurements are valid, we con- 
clude that at energies e ~ 20 GeV the suppression factor 



of the positrons is / 



3, consistent with the results 



of the e + measurements presented in section § [3j 

5. COMPARISON TO DETAILED MODELS 

In this section we describe common models that are 
used in the literature. In § 15.11 we describe the leaky 
box and disk-halo diffusion models and write down sim- 
ple expressions for the grammage in these models. We 
show that detailed models, with parameters chosen to 
fit compositional data, have a grammage that agrees 
with Eq.(fT0|) at high energies. This is done by an- 
alyzing the r e sults of the parameter scan made given 
IMaurin et all (|200lD . In § HTJ we write down expres- 
sions for the loss suppression factors. We show that our 
analytic estimates, using the model's suppression factors, 
reproduce the results of the detailed calculations o f the 
e + flux carried out by IMoskalenko fc Strong! (|1998l ) and 
IDelahave et al.l (120081) to within a factor of ~ 2. As 



commonly claimed ( Adriani et al.1 [2009a). these results 
contradict PAMELA measurements. We highlight the 
assumptions built into these models, that lead to this 
disagreement, and argue that they are not supported by 
observations. 

More details of specific models can be found in 
§ [ A] and in the extensive literature (for reviews 
seeTGinzburg fc Ptuskin] 119761: IBerezinskii et al.l 119901: 
IStrong et al.ll2007D . 

5.1. Grammage in common models 

The CR grammage defined by equations (dJ-® is use- 
ful as long as the conditions for the validity of these equa- 
tions are fulfilled. As these equations are established 
observationally, they must be approximately satisfied by 
any model that is consistent with compositional CR mea- 
surements. These equations will not be valid for a general 
choice of model parameters, in which case the grammage 
is not well defined. For concreteness, we define the gram- 
mage in a general model using Eq. <j5j) , for particles with 
no losses, even if the model does not satisfy equations 
([S])-© exactly. When necessary we specify the distribu- 
tion of sources assumed for this purpose. 

In our calculations, we neglect changes in the particle 
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energy resulting from assumed reacceleration or convec- 
tion. Radiative energy losses of positrons are treated in 



5.1.1. Leaky Box Model 

In this model, the ISM and CRs are distributed ho- 
mogenously in the g alactic disc ('box'), from which 
CRs constantly 'leak' (jGinzburg fc PtuskinlH976l ). This 
model has two free parameters: an effective averaged ISM 
density plbm and a rigidity dependent CR escape time 
T csc (e/Z). The steady state equation for stable particles 
that do not suffer spallation is 



Equation (|35l is equivalent to equation ([5]), with 



X n 



PLBM T csc c - 



(35) 



(36) 



As this model has a uniform composition for all particles, 
primary and secondary, the effect of spallation is exactly 
described by equations ((7j)-<(9j) . 

5.1.2. Diffusion models 

It is commonly assumed that the propagation of the 
CRs is described by a diffusion equation. Models with 
varying assumptions regarding the spatial distribution 
of the diffusion coefficient, sources and boundary condi- 
tions of the propagation volume were studied. In many 
models, the CRs diffuse in a CR halo that surrounds the 
Galaxy, typically a cylinder of radius R > 15 kpc and 
semi-height L ~ few kpc, with a rigidity dependent dif- 
fusion coefficient, while the ISM gas and the CR sources 
are distributed in the galactic gaseous disc with semi- 
height h ~ lOOpc. The CR density for all CR species is 
taken to vanish on the boundaries of the halo. Typically, 
the diffusion coefficient is chosen to be spatially homoge- 
nous and depends on rigidity as D = Do(e/10ZGeV) , 
with D ~ 10 28 cto 2 s _1 and S w 0.5. 

In many examples it is assumed that R ^> L ^> h, in 
which case the radial boundary can be taken as infin- 
ity and the disc as infinitely thin. In this case, we are 
left with a one dimensional model where the distribu- 
tion of source functions of all CRs is the same, namely a 
delta function in z. In this limit, which serves as a good 
approximation for most purposes, the conditions ^) arc 
exactly satisfied. The density of a particle i which suffers 
no losses at a distance z from the disc is given by 



n{z,e) = |^(1 - z/L), 



(37) 



where qi is the generation rate per unit disc area. Equa- 
tion (1371 is equivalent to equation ([5]) when applied to 
the local position z = with the CR grammage given by 



A CS c — X d [ sc Lc/(2D), 



(38) 



where X^isc ~ 10 _3 g cm -3 is the grammage of the 
gaseous disc and the generation rate per ISM mass is 
given by Qi/pisM = <Zi/Adi sc - The effect of spallation 
can be accurately taken into account by equations ([7])- 

©■ 

To illustrate the robustness of eq. ©, as well as to 
demonstrate the calculation of the CR grammage in a 
more complicated set-up, we next consider the detailed 



study of diffusion models made bv lMaurin et al.l (|2001h . 
The results of this study are widely used in recent pub- 
lications. In this study, a x 2 analysis over B/C data 
was made for a wide range of values of the different pa- 
rameters of the model. The model assumes a cylindri- 
cal geometry with radius R — 20 kpc, halo semi-height 
1 kpc < L < 15 kpc, a diffusion coefficient given by 



D(e) = K f3(pc/Z GeV) s 



(39) 



a uniform convective wind with velocity V c , and reaccel- 
eration parameterized by an Alfven velocity Va- 

The parameters Kq,L and S were found to be highly 
correlated. The low x 2 regions were found to be de- 
scribed by the following relation 



K /L = g(L)/f(6), 



(40) 



and the numerical values of f(5) corresponding to several 
values of S were sp ecified. This result ca n be seen in fig- 
ure [7] [taken from (jMaurin et alJl2001Xt )]. The specified 
values of f(8) are shown in figure |8] as black crosses. 

In terms of the analysis presented in the current paper, 
Eq. (j40]) is explained as follows. The reported correla- 
tions simply reflect the fact that, in order to fit the mea- 
sured B/C values, the model's grammage has to agree 
with the measured result (fTUf . More precisely the gram- 
mage constraint must hold at high energies, where the 
energy changing effects of reacceleration and convection, 
that are not taken into account in dZ])-©, are negligible. 
In fact, the correlation (j40j) is determined by the require- 
ment that the value of the grammage at the highest en- 
ergy point in the B/C measurements, e = 35 GeV/nuc or 
pc/Z « e/Z « 75 GeV, be constant and approximately 
equal to the value given by (ITU1) . 

To illustrate this statement, we superimpose onto fig- 
ure [7] contours of constant X esc calculated at fixed rigid- 
ity, e/Z — 75 GeV. The grammage was calculated for 
particles with no losses, for a cylindrical propagation halo 
with an infinitely thin disc. For the disc grammage, we 
used Xdi sc = 200 pcx 1 cm~ 3 x 1.3m p corresponding to 
the 200 pc wide gaseous disc with nn = 0.9 cm -3 and 
"■He = 0.1 cm~ 3 . The halo geometry, R = 20 kpc with a 
solar location at r s — 8 kpc, an d the diffusion coeffi cient 
of Eq. (|39j) were the same as in (jMaurin et al.ll200l . For 
simplicity, the source distribution was taken constant in 
the disc. The grammage for such a configuration, with a 
given value of r s /R, can be written as 



X csc = X disc Lc/{2D)g{L/R) cx e 



-6 



(41) 



where g(L/R) parameterizes the deviation fro m the one 
dimensional model and is given by [using 



](L/R)= 2 -§Y, J o{^ r i 

k=l 



tanh (vkjj) 



(42) 



where J m are the bessel functions of the first kind of 
order m and are the zeros of Jo- Eq. (|41[) can be 
written as 

c^disc g{L/R) 



K /L = 



2X CSC (e/ GeV) 5 ' 



(43) 



reproducing Eq. (|40[) with g oc g. As can be seen in 
figure [7] the reported correlation, Eq. (|40"]) . is explained. 
The non trivial behavior of g(L) reflects the fact that 
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Iso-)j contours for B/C (); <40) 




0.002 0.004 0.006 0.008 0.002 0.003 0.004 0.005 

K /L\kpc Myf'\ ffi) xK,/L \kpc Myf'] 



Fig. 7.— Figure 7 from UMaurin et al.H2001Xl l. The filled gray 
areas on the left panel display contour levels of \ < 40 for the indi- 
cated values of the diffusion power-law index 8. The same contour 
plots are shown on the right panel with the values of Ko/L scaled 
by f(8) (see also text and figure [B). Plotted on top of the figure 
are curves with a given grammage at the rigidity e/Z = 75 GeV 
(e = 35 GeV/nuc), corresponding to the highest energy bin in the 
B/C data (see text). The red solid lines, green dashed line and 
blue dash-doted line correspond to X eS c(75 GeV) = 3.45g cm -2 , 
X csc (75 GeV) = 3.2g cm" 2 and X csc (75 GeV) = 3.9g cm" 2 re- 
spectively. 

the radial extent of the cylinder cannot be neglected in 
the calculation for L > 5 kpc. For smaller scale heights, 
the grammage is equal to the expression in Eq. (|38[) . 
The function f(5), normalized to unity at 6 = 0.6, is 
given by f(6) = (e/ GeV) 5 " 06 w 75 5 " ' 6 . This is in 
good agreement with the required values of f{5) as can 
be seen in figure |S] 

As can be seen in figure [3 the allowed parameters 
have a grammage of X csc (e/Z = 75 GeV) = 3.5 x (1 ± 
0.1) g cm~ 2 in agreement with the value expected using 
Eq. (HQ), X csc (e/Z = 75 GeV) w 3.2 g cm" 2 . 

5.2. Decay and energy loss suppression f s (s) in 
common models 

The calculation of f s _ e + (s) for the above models is 
reviewed in § [AJ The results, along with order-of- 
magnitude estimates are given here. It is useful to 
parameterize the num ber density of CR specie i as 
(jBerezinskii et al.lll990D 

rii = qiT e s,i/V e s,i, (44) 

where is an effective volume probed by the CRs, q~i 
is the total production rate in the effective volume (net 
production in case of particles experiencing spallation) 
and T e s t i is an effective production time. In general, en- 
ergy losses and decay affect both T c s^(e), as the losses 
limit the available time for particle accumulation, and 
KfT,i(s), which is limited by the distance covered by par- 
ticles before they decay or lose a significant part of their 
energy. 




0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 

5 



Fig. 8.— The Ko/L scaling function /(<5) UMaurin et all 12001] . 
black pluses, see also text and figure [71 . f(8) is normalized to 1 
at 8 = 0.6. /(<5) = 75* -6 (solid red line) is required in order that 
the grammage at the rigidity e/Z = 75 GeV (e = 35 GeV/nuc), 
corresponding to the highest energy bin in the B/C data remain 
unchanged. 

In the leaky box model, Qi is uniform and thus 
(ji/V e s,i = Qi is not affected by losses. The production 
time scale is r csc cx X csc in the limit of no losses, and 
approaches r = t c (t ^) in the opposite limit of significant 
energy losses (decay). As a result, we have 

/.(e) ~ r( e )/r eBC ( e ). (45) 

For decay, this equation is exact with r = r^. For energy 
losses, the relation is [see Eq. (jA4|) ] 

f s ,e+( £ ) = ^r c ( £ )/r esc ( £ ). (46) 
7-2 

where 7 is the observable spectral index. The full so- 
lution allowing for moderate losses is given in § [A] As 
can be seen, the suppression factor for decaying particles 
with a decay time Td is similar to that of positrons with 
a cooling time r c if r c ~ Td- This provides an illustration 
for the general argument given in § 14.11 

5.2.2. Disc halo diffusion 

In the one dimensional disc-halo model, the effective 
volume should be replaced by an effective scale height 
L c g. In the case of no losses, the effective height is sim- 
ply the halo size, while the effective generation time is 
roughly the escape time r csc defined by the relation 

L - VDt cbc . (47) 

In the case of significant energy losses (decay), the ef- 
fective scale height is the distance that CRs propagate 
before losing a significant part of their energy (decaying) 
and is given by 

L cff ~ VDW, (48) 

while the effective generation time is simply r = T c (jd)- 
Using Eq. (|4"3|). we thus have 



5.2.1. Leaky Box Model 



(49) 
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It is useful to parameterize f s by: 



fs — Cms- 



(50) 



where Cuie is a dimensionless coefficient that depends 
weakly on the injection spectrum and on the functional 
form of the diffusion coefficient. As shown in § [Aj for 
decaying particles Cms = 1. For positrons with a power 
law injection Q(e) oc e~ 7i and a diffusion coefficient that 

oc 



', Cniff is 

(jA19|) . In the range 7, = 2.5 — 3 and 
the values of Cbiff vary between 0.7 and 



depends on ene rgy a s D oc e X c 
given by Eq. 
6 = 0.3- 0.7, 
0.9. 

Using equations (|49)) . (|27 |) and (|45|) . we see that the de- 
pendence of the cooling suppression factor f se + on en- 
ergy is f s ,e+( e ) e^ -1 ^ 2 , for the disc-halo and 
Leaky Box models respectively. 

In figure [9] we compare the results of the e + 
flux, with the loss given by equation (|50|) . to re- 
sults of calculations g i ven by iDelahave et al.l (|2008() and 
iMoskalenko fe Strong] (|1998h . We used equations (1201) , 



(I50|) with the e + generation cross sections, the pro- 
ton spectra, the disc height and the diffusion coefficient 
adopted from these r eferences. For the comparison to 
(jDelahave et al.l |2008) we calculated the grammage us- 
ing Eq. (|38]>. with A disc = 200 pcx 1 cm" 3 x 1.3m„. 
For co mparison with the results oflMoskalenk o fe Strong! 
(p9l we used A csc = 8.7(e/10 GeVZ)" 6 normal- 
ized to be equal to the value of Eq. (flO|) at e/Z = 
10 GeV and with the pow er law index taken from 
(IMoskalenko fc Strong! fl998l) . 

As can be s een in figure El our ca lculation agrees well 
with that of (IDelahave et al.ll2009f) and up to a factor 
of - 2 with that of (jMoskalenko fe Strong] |l9 98| ) . We 
note that the good agr eeme nt between the calculations o f 
IDelahave etail (|2009h and IMoskalenko fc Strong] (119981 ) 
is necessarily somewhat accidenta l, as the cross sections 
for pr oducing positrons used by IMoskalenko fc Strora 
(119981) are larger than those used by (Dclahav e et al 
2009) by about 40%. 

We next compare the results of diffusion models with 
the PAMELA measurements using the e + + e" measure- 
ments as done in § 13.3! We assume a suppression factor 
/s,e+ = 0.33(e/10 GeV) - ' 25 with the index correspond- 
ing to 6 = 0.5 and normalized to agree with the e + mea- 
surements and the estimate based on radioactive nuclei 
at e ~ 20 GeV. The expected e + flux, divided by e + + e~ 
measurements is compared to the results of PAMELA in 
figure HUl 

As can be seen, the expected flux based on this model 
disagrees with the PAMELA result. We emphasize that 
this inference does not depend on the precise amplitude 
of the e + production cross sections, or the exact value of 
6. This implies that if PAMELA is correct, these simple 
diffusion models are ruled out. 

The fact that the suppression factor is decreasing with 
energy is based on the fact that, in these models, the 
escape time r esc oc e~ s decreases with energy slower than 
the cooling time r c oc e -1 . To see this, note that the 
escape time is related to the grammage in these models 
by 
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Fig. 9. — The interstellar e + flux in diffusion models. The filled 
region is the expected e + flux in case there are no losses as in 
fig [3] The soli d green and dashed blue lines are the results of 
the models in ||Dclaliavc et al. 2008) [ with t he MED propagation 



parameters] and (Moskalenko & Strong 1998) [as parameterized by 
(Baltz & Edsio 1999)] respectively. 
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Fig. 10. — The positron fraction as measured by PAMELA com- 
pared to the expectation for diffusion models. The expected e + 
flux in the model is divided by the e + + e~ measurements. The 
loss suppression is taken to be / e + = 0.33(e/10 GeV) - 25 corre- 
sponding to diffusion with D oc e 0,5 [see J 50 1 ) 1 and normalized to 
agree with the measurements at 10 GeV. The observational data 
is the same as in fig [4] 



V csc ~ V c li sc CT e sc/^- / ; 



(51) 



where we used Eqs. (13"51) and (jr?) . This has a simple 
interpretation - the distance a particle travels before es- 
caping divided by the confining scale is similar to the 
number of disk crossings. Using (|51l) . the fact that the 
grammage X esc changes as e~ s and that L is assumed to 
be constant, we see that r osc changes as r csc oc e~ s . 

A crucial assumption here, that may be wrong, is that 
L does not depend on energy (e.g. lDogiel et al.|[l993l ). In 
fact, a scale height decreasing with energy as L{e) oc e~ Sl 
would imply an escape time changing as r csc oc e~ s ~ Sl 
that may be declining faster than r c if 8l > 1 — 6. Note 
that to a good approximation, the estimate for the anti- 
proton and positron fluxes can be obtained in such a case 
by replacing L in equations (|38|) . (|50|) and (|A18|) with 
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L(e). To see this, note that anti-protons and CR nuclei 
do not change their energy during propagation, imply- 
ing that these equations are directly applicable, while 
positrons and radioactive nuclei are not affected by the 
boundary conditions as long as they are strongly sup- 
pressed by losses. The expected positron fluxes are af- 
fected by the energy dependence of L indirectly, through 
the energy dependence of D(e) inferred from the gram- 
mage measurements. 

Alternatively, r c may deviate from the £ _1 dependence 
if the energy losses strongly depend on the position in the 
halo. 

6. CLAIMS THAT A PRIMARY POSITRON 
SOURCE IS NECESSARY 

The claims made in the literature, that the PAMELA 
data requires a primary e + source (e.g. lAdriani et al.l 
I2009al IMorselli fc Moskalenkdl2008l) . are based on one 
of two different lines of reasonings: 

• The electrons are assumed to have the same pro- 
duction spectrum as the protons, and to suffer the 
same energy losses as the positrons f se - = f s . e +- 
Under these assumptions, the e~ spectrum is given 
(up to an energy independent factor) by 

n e - (e) oc f SiC + (e)n p (e), (52) 

while the e + spectrum is given (up to a factor that 
depends weakly on energy) by 

n e +(e) ex f S}( ,+ X csc (e/Z)n p (e), (53) 

resulting in a e + to e~ ratio, 

n e+ (e)/n e - (e) ex X csc (e/Z) ex £~ - 5 , (54) 

that is independent of the unknown energy loss 
term. The observed slightly rising positron frac- 
tion is in contradiction with this expectation and 
thus a primary source is claimed to be necessary. 
This argument does not rely on direct e~ measure- 
ments. 

• The e + flux, including the energy loss suppression, 
is calculated within a specific propagation model. 
The prediction of the model is divided by the mea- 
sured e~ + e + flux and claimed to disagree with the 
positron fraction [as discussed in § 15.21 see fig [TU] . 

The first line of argument relies on two unestablished 
assumptions, namely that the electrons have the same 
injection spectrum as the protons and that electrons and 
positrons suffer the same cooling suppression. The ac- 
celeration of electrons and their escape from the sources 
are far from being understood. Even if CRs are acceler- 
ated in SNRs, as commonly assumed, and the electrons 
are accelerated with a spectrum similar to protons at a 
given shock velocity, the integrated spectrum of electrons 
leaving the SNR may be non trivial, due to energy losses 
and possible different acceleration efficiencies at differ- 
ent shock velocities. In addition, the electrons may suf- 
fer different losses than positrons on their way to earth 
as their sou rces are distributed differently then the e + 
sources (e.g. lShaviv et al.l l2009; Piran ct al. 2009). 



The second line of argument relies on the details of the 
specific propagation model that is being used (see discus- 
sion in §[5]). The declining suppression factor, which is in 
contradiction with the PAMELA measurements, results 
from the fact that the escape time drops with energy 
faster than the cooling time in the considered models. 
Such a relation between the escape and cooling times is 
not based on observations and can be modified in alter- 
native models (see § I5.2[) . 

7. PREVIOUSLY SUGGESTED MODELS WITH A 
SECONDARY ORIGIN FOR THE POSITRONS 

In § [3] and § E] it was shown that the e + and p mea- 
surements are consistent with model independent con- 
straints of a secondary origin for these particles. How- 
ever, when the commonly used assumptions regarding 
the e~ and e + generation and propagation described in 
§ [5] are adopted, a conflict with the PAMELA results 
arises. In this section we discuss previously suggested 
specific models, which do not adopt the assumptions de- 
scribed in § |6l and which were constructed to reproduce 
the PAMELA measurements with positrons being of sec- 
ondary or igin. 

In (iBlasll [2001 IBlasi fc Serpicol [200l 
iMertsch fc Sarkarl 120091 ) it was suggested that the 
secondary positrons are mainly generated inside SNRs 
resulting with a positron fraction that is flatter than 
in the usual picture of secondary CR production. To 
estimate the contribution of SNRs, the amount of 
grammage traversed by the primary particles during 
the lifetime of the SNRs should be compared to the 
grammage traversed in the Galaxy. The grammage 
traversed in the SNRs is approximately given by 

A^SNR ~ 4tgNR,CpiSM 

n no ( *SNR \ "■ISM _ 2 /„i 

~ °- 08 Tni — i =3§ cm ' 55 

ylO 4 yr ) 1 cm 6 

where £snr ~ 10 4 yr is the age of the SNR and 4pisM 
is approximately the post-shock density. This should be 
compared to the grammage traversed by the primary pro- 
tons at e p ~ lOe through the Galaxy given by 

V Gal ~ 3(e p /100 GeV)- - 5 g cm" 2 . (56) 

As can be seen, the contribution of anti-protons and 
positrons generated in SNRs is negligible for 100 GeV < 
e p < 300 GeV, where the Gala ctic grammage i s mea- 
sured. This is illustrated in ()Blasi fc Serpicol 12009ft 
for anti-protons in the corresponding energy range 
10 GeV < e < 30 GeV. In disk+halo models, the rel- 
ative contribution of the SNRs, compared to the galactic 
contributions for positrons and anti-protons is the same, 
regardless of the different cross sections or the losses. In 
fact, the relative contribution is equal to the ratio of the 
SNR secondary production rate to that in the ISM. It is 
therefor not clear how the positrons at e ~ 10 GeV turn 
out to be d ominated by the contribution from SNRs in 
(|Blasill2009D . 

The effect of acceleration of these secondaries inside 
the SNRs, assuming diffusive shock acceleration theory, 
is limited to the secondaries produced during the acceler- 
ation time. For these, the acceleration amplifies the den- 
sity at a given energy by a factor of a few (acceleration 
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of a e~ 2 spectrum results in a e~ 2 spectrum, somewhat 
upshifted). This is an order unity effect, relevant only at 
the highest energies, where the acceleration time is sim- 
ilar to the SNR age, and where the primary flux begins 
to cut off. 

We note that at e p > 1 TeV, corresponding to sec- 
ondary anti-protons and positrons at e > 100 GeV, 
the contribution depends on unmeasured grammage and 
there is more room for an SNR contribution. In addition, 
if SNRs explode in dense, nisM ^> 1 cm~ 3 regions, the 
contribution may be important. In any case, any such ef- 
fect should be observed for positrons and anti-protons at 
roughly the same energy, and for nuclei spallation prod- 
ucts at higher energies. 

In (|Shaviv et al.l 12001 iPiran et al.l [20091 ) it was sug- 
gested that primary electrons are generated in a nonuni- 
form distribution of sources that follow the spiral arms, 
while positrons are generated throughout the gaseous 
disc. This is an interesting example for the effect of a 
non trivial primary source distribution resulting in a dif- 
ference in the effects of energy losses of electrons and 
positrons. However, the e~ flux is constructed in this 
model to agree with the e + + e~ measurements, while 
the e + are produced roughly uniformly in the disc, simi- 
larly to common diffusion models. It is therefore unclear, 
how the positron fraction obtained in this model is dif- 
ferent than that of ordinary diffusion models, which are 
fitted to the same e + + e~ measurements and in which 
the positrons are generated and propagated in roughly 
the same way (see § I5.2[) . 

Note that this model was not confronted with compo- 
sitional CR data. In particular, it is expected that the 
heavy primary element contribution from the spiral arms 
will be exponentially suppressed due to spallation losses 
i n this model, in contrast t o the observed abundances. 

ICowsik fc Burchl (2009a. b) outlined a model in which 
at low energies particles spend most of their energy de- 
pendent time in the vicinity of the sources before es- 
caping into the Galaxy, where the residence time is 
energy independent. In this model, the primary ob- 
served spectra are the same as the source injection spec- 
tra and the accumulated grammage flattens off with 
energy above a threshold energy e t . By choosing ap- 
propriate parameters, the model can account for the 
observed positron fraction while being consistent with 
nuclei spallation products at e < 10 GeV/nuc. This 
model seems to disagree with seconda ry to primary ra- 
tio measurements at higher energies [jE ngclman et al. I 
119901: iBinns et all Il988t lAhn et all 12008ft . see however 
(jZatsepin et all 12009ft ] A break in the p spectrum at 
~ 10 GeV is expected in such a model. Thus it is cru- 
cial to compare the prediction of these models with p 
observations. 

A feature co mmon to the above model and to the model 
suggested by (Blasi 2009j), is that the spectrum of the 
primary CRs in the regions where the secondaries are 
produced is different than the spectrum at the solar sys- 
tem. These are interesting examples where deviations 
from the expressions for the e + and p fluxes ([TH]) and 
(|20|) are expected, due to the fact that X csc is deter- 
mined by observations of spallation products which are 
generated at rigidities similar to their primaries, while 
the positrons and anti-protons are generated at lower 
energies. For these models, the appropriate grammage 



to consider is the grammage evaluated at the primary 
energy. Note, that the ratio of e+ to p fluxes Eq. ([23]). 
is insensitive to such spectral variations of the primary 
CRs, as illustrated in figure El 

8. SUMMARY AND DISCUSSION 

In this paper we presented a critical analysis of the 
observational constraints on, and the theoretical mod- 
eling of, aspects of CR generation and propagation in 
the Galaxy, which are relevant for the interpretation of 
recent e + and p measurements. 

In §[2] we reviewed the main observationally established 
properties of CR propagation. We focused on the anal- 
ysis of stable CR nuclei measurements and argued that 
they can be accurately described by simple equations, 
(O-© that are valid for a broad range of models, namely 
models which satisfy the conditions ([2]). The basic princi- 
ple behind these equations is that the CR density should 
be proportional to the CR generation rate as long as the 
CR composition is uniform. The main challenge in the 
analysis of secondary CRs, including anti-protons and 
positrons, comes from the fact that CRs suffer during 
their propagation significant losses, which are different 
for different species. Nuclei suffer spallation and decay, 
while positrons suffer energy losses due to Inverse Comp- 
ton and Synchrotron emission. The spallation losses of 
stable nuclei can be absorbed into a net production rate, 
see Eq. ©. 

The CR grammage A csc , see Eq. (fit)]) , extracted from 
the stable CR nuclei measurements, allowed us to write 
down in § |3] analytic, model indepe ndent constraints for 
the expected secondary p flux (a la iGaisser fe Schaefer I 
Il992t iSimon et all 119981 equivalent to Leaky Box and 
thin-disc+halo models), and for the secondary e + flux, 
including an accurate estimate of the p flux, Eq. (|18p. 
and an upper limit for the e + flux, Eq. (|23j) . obtained 
by neglecting the e + radiative losses. We have shown 
(§ [5]) that the simple analytic expressions reproduce the 
results of detailed calculations, based on popular mod- 
els for CR propagation, including t he leaky box model 
and disc+halo diffus i on models (e.g. iMaurin et alj[2"001b 
iDelahave et~aDl2008l : iMoskalenko fc Strong! 1 19981 ). Anti- 
proton and e + measurements available up to e ~ 100 
GeV, including the PAMELA positron fraction measure- 
ment, were shown to be consistent with our analytic es- 
timates in figures [2]l5] 

The claims that the measured PAMELA positron frac- 
tion requires primary e + sources were shown in § [6] to 
be based on assumptions, that are not supported by ob- 
servations. If PAMELA results are correct, the fact that 
they are in disagreement with the results of models based 
on these assumptions simply imply that the assumptions 
are not valid. In particular, the models which are ruled 
out include (a) models in which the primary electrons are 
injected with a source spectrum equal to that of protons 
and in which the electrons and positrons suffer identi- 
cal losses, and (b) thin-disc+halo diffusion models, with 
an isotropic diffusion equation with energy independent 
halo size. 

Positron flux measurements may be used to infer the 
suppression factor f s e + (e) [see Eq. (JTJ] due to radiative 
energy losses. Existing e~ and e + data were shown in §[3] 
to imply that f se +(e) ~ 1/few in the PAMELA energy 
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range, and that f Sye +(s) is indicated to be slightly grow- 
ing with energy, see figure [4] and [5] and the discussion 
following them. In § [4] we showed that measurements of 
radioactive unstable secondaries allow an indirect esti- 
mate of / S)e + ~ 1/3 at a single energy, ~ 20 GeV (§ I4.2|I . 
consistent with the e + measurements. The e + measure- 
ments at higher energies, e > 20 GeV, provide new infor- 
mation on the temporal behavior of CR propagation that 
is not currently accessible by other types of experiments. 
Measurements of the e + to p flux ratio are particularly 
useful for this purpose, as well as for challenging sec- 
ondary production models, as this ratio is insensitive to 
model details (see figure [6]). 

In section §[7]we discussed specific models, in which the 
positrons are of secondary origin, that were s uggested to 
explain recent e + abundance measurement s (|Blasill20Q{l 
IShaviv et al l [20091 : ICowslk fc Burchll2009aD . We showed 



that these models appear to be in conflict with con- 
straints derived from CR measurements. It is important 
to emphasize, that there is a variety of models that can 
be constructed, in which both e + and anti-protons are of 
secondary origin, and which are consistent with the con- 
straints arising from CR observations. The main point 
that should be addressed in such models is the energy 
dependence of the e + suppression factor f se + at ener- 
gies exceeding e > 20 GeV, which is indicated by the 
PAMELA measurements to be slightly rising with en- 
ergy. One example, that could result in such a behavior, 
is to have energy depen dent boundary con ditions to the 
transport equation (e.g. iDogiel et al.lll993l see 



We thank W. R. Webber, T. Volansky and P. Salati 
for useful communications. This research was supported 
by ISF, AEC & Minerva grants. 



APPENDIX 

CR PROPAGATION MODELS 

Below we provide the calculations quoted in jj [5] Calculations along similar lines can be found in the literature (e.g. 
iGinzburg fc Ptuskin|[l97llBerezinskii et al.| [l990t Nevertheless, we find it convenient to attach in full the treatment 
of specific problems encountered in the paper. 

Cooling suppression of positrons in the Leaky Box Model (LBM) 
For positrons, the LBM equation reads 



Q e + = + d e (in e +) 

To.a,r. 



where — e — e/r c oc e 2 is the energy loss rate. Solving Eq. ()Alj) we obtain the suppression factor due to cooling, 



r c (e) JIT d£ 'Qe+( £ ') e *P 


f£ de" r c {s") 
Je £"r c3c (e") 


Tcsc(e) £<9e+ 


(e) 



In the limit t c <C t csc , the exponential term can be omitted and we get 



fs,e+ (e) 



T esc (e) Q e +(e) 



7i - 1 t C! 



(Al) 



(A2) 



(A3) 



where the last equality holds for a power law injection spectrum, Q e + (s) oc e 7i . In the limit r c <C t csc , the injection 
index 7, is related to the observed index 7 by 7, = 7 — 1 and so Eq. (|A3|) can be written as: 



f, 



1 



7 - 2 t c . 



(A4) 



Diffusion model with Halo and thin disc 

Stable secondaries with no energy losses, cylindrical model — For the case considered in §0 it is useful to perform a bessel 
decomposition in the radial coordinate while solving for the z coordinate directly. The solution for the CR density is 
given by 

tanh [vk^jA cosh {vkj^J - sinh {^kj^j 



n{z, r; e) 



R T ( T \ 1k{e) 



D(e) 



k=l 



RJ 2v k 



(A5) 



where qk are the coefficients in the bessel decomposition of the source surface density function q(r;e), defined on the 
gaseous disc. We assume a source function which is uniform across the disk q(r; e)—q(e) 1 which has the following 
decomposition 

= ~^-y (A6) 

The grammage is defined by the equation n = X esc Q / (pis m c ) [Eq. ©]. Here the local generation rate per unit mass 
is given by Q/pism = q{r s ; e)/A^di SC - Using (|A5|) and (|A6|) . the local density can be expressed as 



Die) 



k = l 



tanh (vkji) 

R) vpiiyk) 



(A7) 
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The grammage for this model is therefore 



Lc 2R 
2D X 



tanh \ Vk 



k=l 



Rl v\Jx{y k ) 



(A8) 



In the limit L<R the summed factors tanh(^L/i?) — > VkLjR, and the result Eq. (l38l) of the one dimensional model 
is reproduced. 

Decaying particles, one dimensional limit — For decaying particles with observer frame decay time r(e), the solution to 
the diffusion equation is 



tanh 



'Dt 



cosh 



where q, D, r depend on e. The suppression factor defined at z = is 



fs 



'Dt 



tanh 



L 



— sinh 



'Dt 



'Dt 



Dt J ^Tdt^l L 



(A9) 



(A10) 



Continuous energy losses, one dimensional limit — We next consider particles that suffer energy dependent cooling e[e). 

It is useful to express the steady state density in terms of the time dependent number density of particles n(t, e)dt 
that would result if the source generated particles at all energies and throughout the disk at an interval dt at time 
t = 0by 



n(e) 



n(t, e)dt. 



It is straight forward to show that 



n(t,e) = q(e s )—n [t(t 1 e s )] 



where es is the energy a particle had at time t = in order to be at energy e at a fictitious time t defined by 



t(t,e s )= / D[s(e s ,t')]dt' = T c (e)D(e) 



es/e 



D(£e) i(e) 



(All) 
(A12) 

(A13) 



D(e) ifc) 

and no(t) is the value at z = of the solution to the spatial diffusion equation dfno = d^no, with initial distribution 
no(t = 0,z) = 5(z) and a constant diffusion coefficient equal to unity. Using Eq. (lAlljl . the steady state distribution 
is given by 



n(e) = q(e)T c (e) 



l 9(e) 



■n (t)d£. 



(A14) 



For source and diffusion coefficient with a power law dependence on energy, q(e) oc e li and D(e) oc e s , and taking 
e oc e -2 we find 



t{ss) = r c (e)D(e) 



£, s ~ 2 d£_ 



1 



[1 - (es/e^T^Die) 



1-5 



and 



n = qr c 



1 - S 



(i-e- L )T C D\dt 



where all physical quantities (n,q, D,t c ) are to be evaluated at e. In the limit of strong suppression we find 



qr c 



vr~5 r 7< (i-^ _1 )" 1/2 ^ 



\J AttDt c 

in agreement with the result in (jDelahave et al.ir2008l ). The implied suppression factor is 



fs 

where Coiff is a dimcnsionlcss coefficient given by 



CDifF" 



L ' 



(A15) 
(A16) 

(A17) 

(A18) 

(A19) 



In the range ji = 2.5 — 3 and 6 = 0.3 — 0.7 the values of this function are between 0.7 and 0.9. 
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